Fluctuations in a superconducting quantum critical point of multi-band metals 
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In multi-band metals quasi-particles arising from different atomic orbitals coexist at a common 
Fermi surface. Superconductivity in these materials may appear due to interactions within a band 
(intra-band) or among the distinct metallic bands (inter-band). Here we consider the suppression 
of superconductivity in the intra-band case due to hybridization. The fluctuations at the supercon- 
ducting quantum critical point (SQCP) are obtained calculating the response of the system to a 
fictitious space and time dependent field, which couples to the superconducting order parameter. 
The appearance of superconductivity is related to the divergence of a generalized susceptibility. For 
a single band superconductor this coincides with the Thouless criterion. For fixed chemical potential 
and large hybridization, the superconducting state has many features in common with breached pair 
superconductivity with unpaired electrons at the Fermi surface. The T=0 phase transition from the 
superconductor to the normal state is in the universality class of the density-driven Bose-Einstein 
condensation. For fixed number of particles and in the strong coupling limit, the system still has an 
instability to the normal sate with increasing hybridization. 
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INTRODUCTION 



In heavy fermion materials superconductivity can be 
suppressed in different ways. Most commonly, this is ac- 
complished by an external magnetic field, applied pres- 
sure or doping [J-Q, but the point in the phase diagram 
where the critical temperature T c vanishes as a function 
of the external parameters is not necessarily associated 
with a SQCP. For example, in the case of superconductiv- 
ity induced by antiferromagnetic fluctuations close to an 
antiferromagnetic quantum critical point (AFQCP)|9|, as 
the system moves away from the AFQCP, these fluctua- 
tions change from attractive to repulsive and supercon- 
ductivity just fades away [10|. 

Recently, we proposed a new mechanism for destroying 
superconductivity which is specific for multi-band mate- 
rials as heavy fermions [l], [ll|. It occurs because the 
narrow band of quasi-particles which form the Cooper 
pairs is hybridized with a band of normal electrons. As 
hybridization increases by pressure or doping it turns out 
that a minimum attractive interaction is required to pro- 
duce a superconducting ground state. This mechanism of 
suppression of superconductivity has been demonstrated 
using a mean-field approach. Starting in the super- 
conducting phase as hybridization is turned on, T c was 
shown to vanish at a critical value of the mixing (V c ). 
However, the mean-field approach does not include fluc- 
tuations. Although the normal state for V > V c is metal- 
lic, there are no critical fluctuations associated with the 
onset of superconductivity in this approximation. Here 
we go beyond the mean-field approximation to include 
fluctuations. We start in the normal phase, where either 
V > V c or the attractive interaction is not sufficiently 
strong to produce a superconducting ground state. As 
hybridization decreases, or U increases, the system has 



an instability to a superconducting ground state. This is 
shown using a new approach [l2[ based on a perturba- 
tion theory for retarded and advanced Green's functions 
[l3j |. For the single band case our results reduce to the 
well known Thouless criterion 1J] for superconductiv- 
ity. We relate the appearance of superconductivity in 
the multi-band metal to the divergence of a generalized 
susceptibility, very much like the Stoner criterion for the 
appearance of ferromagnetism. 

We calculate the response of the system to a frequency 
and wave- vector dependent fictitious external field which 
couples to the superconducting order parameter. For 
simplicity we consider here the case of s-wave super- 
conductivity. Using an RPA approximation, equiva- 
lent to summing an infinite series of bubble diagrams, 
we obtain a generalized susceptibility x(q,w). At zero 
temperature, the static and homogeneous part of the 
susceptibility diverges at exactly the critical mean-field 
value of hybridization, which destroys superconductiv- 
ity. This divergence implies that even in zero field the 
system can have a finite superconducting order param- 
eter. The condition for the superconducting instability 
can be expressed in the form of a Stoner-like criterion, 
1 — U Xq(V) = 0- I n the case with a fixed chemical poten- 
tial this determines either a critical value of the attrac- 
tive interaction U C (V) above which the system is super- 
conductor or a critical hybridization V C (U) below which 
superconductivity sets in. 

The approach developed here allows to obtain the na- 
ture of the fluctuations close to the SQCP. We find that 
at U C {V) (or V C {U)) and for low frequencies and small 
wave- vectors, the generalized frequency and wave- vector 
susceptibility has poles at real frequencies u> = T)q 2 . It 
turns out that the SQCP is in the universality class of 
the Bose-Einstein condensation 15]. From the nature 
of the quantum critical fluctuations, we obtain the dy- 
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namic critical exponent z using Hertz approach to quan- 
tum phase transitions [16[. This allows us to obtain the 
thermodynamic properties when the system is close to 
the SQCP [17J. Our theory is equivalent to a quantum 
Gaussian approach and since the dynamic exponent turns 
out to be z = 2, it yields the correct description of the 
SQCP for dimensions d > 2. In our treatment the sys- 
tem is metallic in the non-superconducting side of the 
phase diagram. This arises due to quasi-particles which 
remain unpaired close to the Fermi surface. This is dif- 
ferent from the predictions of disordered induced SQCP 
(ljj . where the normal state is insulating at T — due 
to the presence of a gap. In the case of a single band and 
at zero temperature, the Thouless criterion for BCS su- 
perconductivity implies a superconducting ground state 
for any finite value of U. 

It is interesting at this point to comment on the dif- 
ferences between the intra and inter-band problems. We 
have learned from the mean-field approach [ill ] that in 
the former case, superconductivity is destroyed contin- 
uously through a SQCP as hybridization is increased. 
On the other hand for inter-band superconductivity, as 
hybridization increases superconductivity is suppressed 
discontinuously at zero temperature. This occurs at a 
first order quantum phase transition (ill Il9j ] with phase 
separation and metastable regions. For intra-band at- 
tractive interactions as one approaches the superconduct- 
ing state from the normal phase, the instability of the 
normal metal is for a q — BCS state and occurs at 
the same value of hybridization for which the homoge- 
neous BCS state is destroyed. For the inter-band prob- 
lem, the abrupt disappearance of superconductivity sug- 
gests to look for other forms of superconductivity which 
can remain stable in the presence of a large Fermi wave- 
vectors mismatch. For inter-band attractive interactions, 
the normal metal becomes unstable to an inhomogeneous 
q 7^ superconducting state for a value of the Fermi 
wave- vectors mismatch larger than that for which super- 
conductivity is destroyed through a first order transition 
Q. In a previous work, using the same approach pro- 
posed here, we could determine the universality class of 
the quantum phase transition from the normal to a pair 
density wave phase as the mismatch of the bands with 
inter-band interaction increases (l2l | . 

The study of the crossover from weak coupling to 
strong coupling has attracted considerable attention after 
the discovery of the high-temperature superconductors. 
More recently, new experiments in cold atom systems 
have further raised the interest in this problem [20j|. In- 



troducing the concept of scattering length, l/kpa s , [21 1 
we are able to analyze the case with fixed number of par- 
ticles and the crossover between the weak and strong cou- 
pling for the two-band system with hybridization. This is 
necessary since the values of U for which there is a critical 
value of hybridization are in the strong coupling limit. In 
this case the gap equation must be solved simultaneously 



with the number equation and a trivia (ji, A, l/kpa s ) can 
be obtained numerically. We find that the critical value 
of the hybridization increases with the increase in the 
coupling strength tending to a saturation value. 



THE MODEL 

We start with the following Hamiltonian describing a 
two-band system 22j , with a hybridization between them 



and an attractive interaction in one of the bands 
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+ £ WU*r + aj>*0> (!) 

where <zj ff , dj a , b\ a , bj a create and destroy electrons in the 
narrow a-band and the wide &-band, respectively. The at- 
tractive interaction U whose origin we do not specify acts 
only between the electrons in the narrow a-band. The hy- 
bridization V mixes the electrons of different bands and 
can be controlled by external pressure. In order to cal- 
culate the superconducting response of this system, we 
proceed with the approach proposed in Ref. [l2j . intro- 
ducing a wave- vector and frequency dependent fictitious 
field that couples to the superconducting order parame- 
ter 
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( 8 ki + fl .t a 'i). ( 2 ) 



where the frequency u>o has a small positive imaginary 
part to guarantee the adiabatic switching on of the field. 
The response of the system to the fictitious field, will be 
obtained using perturbation theory for the retarded and 
advanced Green's functions [l3[. We start in the normal 
phase where the superconducting order parameter is zero 
in the absence of the external field g. We split the Green's 
functions, normal and anomalous in two contributions. 
The first of order zero and the second of first order in the 
field g. For the anomalous Green's functions, we write, 



« al\a] a , » -> « 4>j ff , » (0) + « 4>j v » 

In the normal phase and in the absence of the fictitious 
field, the relevant zero order Green's functions can be 
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easily calculated [ll|. They are given by, 
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The first order propagators should be obtained more 
carefully due to the time dependence of the external field. 
Let us consider the equation of motion for the first order 
anomalous Green's function in the normal phase, 



«L(*)I«L(*') = -£i*3« at(*)l«L(*') » 
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where a — ±1 para <r =t>4- 

We decouple the higher order Green's function in the 
spirit of the random phase approximation (RPA) to ob- 
tain 
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*Ji« «lMhUO » = -£i««« aL(*)hU*') » 

+(7<5A 1CT (£)« en-a^la],,^) > (0) 
-E ; ^«bL(t)|a] CT ,(0» (1) 
+p5e- lq - r 'e^ 0t « Oi-^CtJIaJ,,,^) » (0) , 
where in the adiabatic approximation, 

6A ia (t) =SA a a a e~ lcl - ri e lUc ' t . (4) 

Next, we perform a Fourier transformation only in the 
time variable t, using that a\ a (t) = J duia icr (ui)e~ luf and 

the same for b\ a (t). Notice that the zero order propaga- 
tors are time translation invariant and depend only in the 
time difference (t — t 1 ). The first order propagators how- 
ever are functions of both t and t 1 . We can also Fourier 
transform in space, using that a\ a {t) = Ek a L(*) elk ' rii 
We get 
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Spatial translation invariance is lost due to the spatial 
dependence of the external field. Now, we go on to write 
the equations of motion for the new generated Green's 
functions. Proceeding like above we arrive at the follow- 
ing system of equations, 

(" + e£)«4>t,_ CT »^ = WAr<<a (fc _ g) _>t,_ CT >>^ o 
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( W + 4)« &1>U >i x) = -v fc < a Ll a l'_ CT >i 1) , (7) 



(w - 6 b fc )« 6 fc - CT |4,_ CT »^ = 14« a fc _ CT | a t,_ CT »^. (8) 

We notice from Eq. 0that inter-band pairing can be in- 
duced by hybridization even in the absence of inter-band 
attractive interactions. In the case of purely inter-band 
interactions, the normal-BCS superconductor transition 
is discontinuous and phase separation occurs 25, l26j |. 
Here, however, these inter-band anomalous correlations 
induced by hybridization do not affect the nature of the 
transition 

The system of equations [MS can be easily solved. In 
particular we get for the anomalous Green's function, 
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Notice that the index (t') becomes redundant since the 
zero order Green's functions are time translation invari- 
ants. The anomalous first order propagator can finally 
be written as, 

«a^\a^ - (1) - 1 v 



(u+u -e k _ ) 



[USA aa + ga] . (10) 



Then, the anomalous first order Green's function can 
be completely determined in terms of the zero order equi- 
librium Green's functions obtained previously. This al- 
lows to use the fluctuation-dissipation theorem [l3j to ob- 
tain the correlation function <5 A aa from the above anoma- 
lous Green's functions, 
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where F w {G(cj)} = -/ dojf(u})[G(oj + ie)-G{u>-ie)] is 
the statistical average of the discontinuity of the Green's 
functions G(uj) on the real axis [Hj]. The function f(u>) is 
the Fermi-Dirac distribution. Defining the susceptibility 
(see Eq|3]) as, 

X a a (q, u ) = -2irJ2 F^{G aa (-"i k)Q aa {io + uj ,k- q)}, 

k 

we can write the response to the fictitious field in the 
form, 
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Then, the change of the order parameter due to the 
fictitious field is given by a generalized susceptibility 
X aa {q,uo), in the form 5A aa = X aa (<L w o)<?- As param- 
eters of the multi-band system in the normal metallic 
phase, such as, the hybridization or the strength of the 
attractive interaction change, eventually the condition 
1 — Uxo a (q,^o = 0) = is attained. This signals an in- 
stability of the system to a superconducting phase since, 
even in the absence of the field g, the order parameter 
may be finite. So, the appearance of superconductivity 
is related to the divergence of a generalized susceptibil- 
ity. If the above criterion is first satisfied for q = 0, 
the instability is towards a homogeneous superconduct- 
ing phase, if it occurs first for a finite q = qa, an in- 



homogeneous superconducting FFLO-like phase 27| sets 
in. This is the case in general for inter-band interactions 
in the presence of a mismatch of the Fermi-wave-vectors 
of the non- hybridized bands as is shown in Ref. [Hj]. 
For intra-band interactions, the instability is towards a 
homogeneous superconducting phase when the condition 
1 - Uxo a {q = 0, w = 0) = is fulfilled [ll|. 



The generalized susceptibility 

Making a change of variables, and writing in a more 
symmetric form, the frequency and wave-vector depen- 
dent susceptibility of Eq. [TT]is obtained as, 
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where. 
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are the energies of the hybrid bands. These bands are 
schematized in Fig. [1] 
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FIG. 1: (Color online) Dispersion relations for nomothetic 
bands. The dashed lines represent the original unhybridized 
bands and the full lines the hybrid bands. For V > Vbcs (see 
text), the region delimited by the cut-off vector kc around the 
Fermi wave- vector Izf of the original bands does not include 
the Fermi wave- vectors Izfa and fcps of the hybrid bands, as 
shown in the figure. 



In the case hybridization vanishes, V — > 0, we get, 

m +q )-f(4- q ) 



hm X g a (2g, Wo ) 



E 



-k+q 



(15) 



and the condition 1 — Uxo a (q — 0, uiq — 0) = 0, yields, 
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which is the BCS gap equation [28|, where V is the vol- 
ume and P = 1/ksT. At zero temperature this Thouless 
condition is satisfied for any finite U due to the logarith- 
mic divergence of the integral. On the other hand, in 
the presence of hybridization V, for q — > e uio — > the 
susceptibility is given by, 



xS a (o,o) =E» 



[i-2/K)] 



/ 62 a 



, b2 1,2, 



(16) 



In this case the condition 1 — U Xo a (q — 0, oj = 0) = 
at T = yields a critical value for hybridization V c be- 
low which the system becomes superconductor. Alterna- 
tively, for a given value of V, there is a critical value of 
the attractive interaction above which the system is su- 
perconductor. This condition turns out to be the same of 
that found previously using a mean-field approximation. 
In this case, however, the system starts in the supercon- 
ducting phase and as hybridization reaches the critical 
value V c it enters the normal phase continuously at a 
SQCP. Further understanding of this instability is gained 
by examining Fig. [TJ It shows the dispersion relations 
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for the original bands (V = 0) and when the hybridiza- 
tion is turned on. For simplicity we consider nomothetic 
bands, e b k — ae^, (a > 1), with e£ = k 2 /2m a — fi (H = 1). 
The integration in /c-space in the generalized susceptibil- 
ity is done in a region of width 2kc around the original 
Fermi wavevector kp (k F /2m a — (J-)- As V increases 
there is a reorganization of the electronic structure. The 
bands repel and the Fermi wave-vectors of the new hy- 
brid bands eventually move out of the region of integra- 
tion (see Fig. Q} avoiding the weak coupling logarithmic 
singularity. The interest of the model is that, even in 
the situation of Fig. [TJ the system can sustain a su- 
perconducting ground state if the attractive interaction 
is sufficiently strong. In these conditions, the equation 
1 — Uxo a (q = 0,o;o = 0) = can be solved numerically 
for finite temperatures and a critical line as a function 
of hybridization, T c x V, ending at a SQCP can be ob- 
tained . Close to Vc the critical temperature vanishes 
like, Tc oc \J~Vc ~~ V with a mean-field shift exponent 
^ = 1/2. 

Solving the equation 1 — Uxc> a (q = 0,w = 0) = 
for T = 0, we obtain the zero temperature phase dia- 
gram of Fig. [2j It shows the critical value of the in- 
teraction versus hybridization. Below a value that we 
called Vbcs the system presents the usual BCS instabil- 
ity and is superconductor for any finite U. For values 
of V > Vbcs, a finite value of the interaction is neces- 
sary for the existence of superconductivity. So, in this 
case, we have a quantum phase transition tuned by the 
strength of the hybridization for a given value of interac- 
tion or tuned by the strength of the interaction for a given 
value of hybridization. The value Vbcs is that for which 
the Fermi wave- vectors of the new hybrid bands coincide 
with the limit of integration around the original Fermi 
wave-vector. For homothetic bands Vbcs — 2~ s favpkc- 
As V — > Vbcs from above, the critical line vanishes as 
the inverse of tanh -1 (V/Vc) (see dot-dashed line in Fig. 

Further knowledge of the nature of the quantum phase 
transitions requires to analyze the frequency and wave- 
vector dependence of the generalized susceptibility, Eq. 
1131 Here we are interested in the metal-superconductor 
transition for V > Vbcs, he., in the case the system en- 
ters the superconducting state by reducing the hybridiza- 
tion or increasing U as along the trajectories shown by 
the dashed lines of Fig. [2] In these conditions, we can 
obtain the real and imaginary part of the dynamic sus- 
ceptibility. The imaginary part turns out to be zero and 
the real part is given by: 



Re[^ a {2q,Lj + ir})\ = 



3 fc c 1 
16tt 3 k F V 
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where F(a) — 10(1 + a), 5kp — jx/avp and k c the cut-off 
vector in k-space. Then, at the critical value of hybridiza- 
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FIG. 2: (Color online) Zero temperature phase diagram for 
the intra-band superconductor for fixed chemical potential. 
In the normal phase, as the hybridization (V) is reduced, for 
finite U, or U is increased for V > Vbcs as indicated by the 
dashed lines, the normal metal becomes unstable to a super- 
conducting phase. This phase transition is in the universality 
class of the zero temperature Bose-Einstein condensation as 
discussed in the text. For V < Vbcs, T = 0, the system is 
a BCS superconductor for any value of U 0. The values 
of the parameters used to obtain the curve are: a = 5 and 
tdc ~ VFkc = 0.001. All units are in terms of the bandwidth 
fi. 



tion the system has propagating bosonic excitations with 
a quadratic dispersion relation given by: 



w 



2m, 



-F(a). 



We can rewrite the dynamic susceptibility as: 
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The critical line V C (U) or U C (V) for V > V BC s is 
shown in Fig. [2] (dot-dashed line). Along this line, the 
dynamic critical exponent associated with the normal- 
superconductor transition is z — 2. This transition is 
in the same universality class of the zero temperature 
Bose-Einstein condensation as a function of the chemical 
potential (l5| . This type of transition occurs in anti- 
ferromagnetic systems at T = where spin-waves can 
also present the phenomenon of Bose-Einstein condensa- 
tion as a function of the magnetic field [15| , It is useful 
to express the results above using the language of the 
renormalization group. The interaction U is a relevant 
parameter at the fixed point, Vbcs, U = 0, and the 
quantum critical behavior all along the critical line in 
Fig. [2] for V > Vbcs is governed by a strong coupling 
fixed point associated with the Bose-Einstein condensa- 
tion of Cooper pairs. Since the effective dimension for 
this transition d e ff = d + z is greater than the upper 
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critical dimension in three spacial dimensions, its critical 
exponents assume Gaussian values. A scaling analysis 
17] using z = 2 and Gaussian exponents allows to de- 



termine the thermodynamical properties of the system in 
particular along the critical trajectory, (V = Vc, T — > 0). 



The critical line in d = 3 approaches the SQCP as [17j . 
T c oc \V - Vb| 2/3 . The contribution of the critical fluc- 
tuations for the specific heat along the critical trajectory 
17j is given by C oc T 3 / 2 . For disordered systems, the 



scattering of electrons by bosons with parabolic disper- 
sion relation yields a resistivity, p oc T and p oc T 3 / 2 in 
d = 2 and d = 3 [2§| , respectively. 



Stability of the superconducting phase 

The superconducting phase we obtain for V > Vbcs 
has interesting features. It has pairs of quasiparticles co- 
existing with gapless excitations at a Fermi surface. In 
some respects it resembles the breached pair supercon- 
ducting phase proposed some time ago by Wilczek and 
Liu [30j. It is then natural to ask about the stability of 
the superconducting phase for large values of hybridiza- 
tion. For this purpose we study the difference between 
the ground state energy of the normal and superconduct- 
ing phases as hybridization increases. The quasi-particle 
excitations in the superconducting phase have been ob- 
tained in Ref. 



Ill] and are given by, 
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where A aa is the superconductor order parameter. The 
free energy due to condensation is given by, 
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ki? —k c 



n h k )dk- 



(A 



aa\2 



u 



(23) 

This vanishes in the normal state and its minimization 
with respect to the order parameter A aa yields the gap 
equation for the intra-band superconductor order param- 
eter. In Fig. [3] we show this energy difference as the sys- 
tem crosses the phase transition line for a fixed V > Vbcs 
with increasing U, as along the vertical dashed line in Fig. 
[21 The phase transition from the normal (A = 0) to the 
superconducting state with increasing U is continuous or 
second order. The quantum phase transitions at the crit- 
ical line U C (V) for V > Vbcs as shown before are in the 



universality of the density-driven quantum Bose-Einstein 
condensation. The superconducting phase is clearly as- 
sociated with a minimum in the condensation energy as 
shown in Figure [3] 

Figure @] shows the energy difference AE as the su- 
perconductor moves from the Bose-Einstein condensation 
regime at V > Vbcs to the BCS regime for V < Vbcs 
with decreasing hybridization. It is clear from the fig- 
ure that this change of regime has a smooth behavior 
and nothing special occurs at Vbcs for a finite value 
of U. As the BCS phase is stable for V < Vbcs, this 
figure with the continuous behavior of the ground state 
along the change of regime implies also the stability of 
the BEC-like phase for V > Vbcs- 




FIG. 3: (Color online) Condensation free energy, Eq. 1231 for 
different values of the interaction. The system with a fixed 
value of hybridization (V > Vbcs) and U/p = 1.1 starts in 
the normal state with a minimum at A aa = 0. As the attrac- 
tive interaction U/p increases, it becomes a superconductor 
through a continuous quantum phase transition which is in 
the universality class of the density-driven Bose-Einstein con- 
densation. This phase transition corresponds to following the 
vertical dashed line in Fig. [2] Values of the parameters and 
units as in Fig. [5] 

Our results are obtained with the superconducting 
pairs being formed close to the original Fermi wave- 
vectors of the unhybridized bands. One may argue that 
as hybridization is turned on, we should integrate not 
around the original kp but around kpB of the new hy- 
brid band in which case mixing would never suppress 
superconductivity. However, in the process of diagonal- 
ization the two-band problem, terms of the type cvw , 

/3 t and ft fft are generated. The operators a\ a, ft , 
P are the creation and annihilation operators of the new 
hybrid bands. Then all types of superconducting corre- 
lations among the hybrid bands appear and there is no 
particular reason to integrate around, say kpB- In order 
to expose the complexity of this problem, we calculate 
the difference between the ground state free energies of 
the superconducting states obtained integrating around 
&f and kFB- This difference is shown in Fig. [5] For the 
BCS type states with V < Vbcs, this is always nega- 
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FIG. 4: (Color online) Condensation free energy, Eq. 1231 as 
the system moves from the BEC state to the BCS state with 
reducing hybridization (horizontal dashed line in Fig. [2]). As 
shown, both phases are smoothly connected. The supercon- 
ductivity minimum just shifts and gets deeper continuously as 
V decreases passing through Vbcs- Values of the parameters 
and units as in Fig. [2] 




formalism in which we are able to reach also the strong 
coupling limit. In this section we study the quantum 
phase transition tuned by hybridization in a range of the 
interaction strength. Previous studies in one-band sys- 
tems have shown that the evolution between these two 



different regimes is continuous [3 lM33j j . For the case of 
a two-band system with d-wave symmentry it was also 
found that this crossover is smooth 34 1 . 

We start with the generalized gap equation to describe 
the weak to strong coupling crossover extended for the 
two-band case with hybridization, jj = Xo a (2g,a>o) In 
this equation Xo a (2<?,u;o) is given by Eq. Q2] with the 
replaced by the dispersion relations in the supercon- 
ducting state, Eqs. |2"0U22"1 

In this case a BCS-type cut-off cannot be used to treat 
the strong coupling limit, in which all the particles must 
interact, and not just those around the Fermi surface 
as in the BCS limit. We thus introduce the concept of 
s-wave scattering length, a s , to regulate the ultraviolet 
divergence in the gap equation that arises due to the 
sum over all energies [21(. It is related to the coupling 
strength U through: 



47ra, 



1 

U 



E 



i 



2(4 



(24) 



and can be positive or negative. When \/kpa s — > — oo 
one obtains the weak coupling regime, while l/kpa s — > 
oo gives the strong coupling limit. 

Eliminating the coupling U from the gap equation us- 
ing the relation above, we find: 



FIG. 5: (Color online) Difference of the energies of the super- 
conducting states integrating in shells of width k c around kF 
and kFB (see Fig[T|, for the same values of V and parameters 
as in Fig(4] Dashed lines are for V > Vbcs and continuous 
lines for V < Vbcs- Notice that the BCS state for V < Vbcs 
has always a lower energy when integrating around the orig- 
inal k F (AE < 0). Also for V > Vbcs, A£ < at the 
minimum of the superconducting state showing the stability 
of the breached pair-like state for U > U C (V) with respect to 
the BCS state obtained integrating around kFB- 



tive showing that integration around kp leads to a lower 
energy. For V > Vbcs, the breached pair type of super- 
conducting state has, at the value of A aa for which the 
energy is a minimum, a lower energy than that of the BCS 
state obtained integrating around fe^g for V < V C (U). 



THE BCS-BEC CROSSOVER 

From Fig. [2] one can see that for values of V slightly 
greater than Vbcs the value of U c is already around Ep, 
so is seems important the study of the transition in a 



47ra< 



1 



2(4 



M) 



(25) 



For the purpose of solving this equation and determine 
V c we can fix the chemical potential or the total number 
of particles. For the second case we need first to obtain 
the chemical potential as function of the new coupling 
parameter \/kpa s . From the propagators in the super- 
conducting phase given in ll| we can find the number 
equation, that guarantees the conservation of the total 
number of particles through the variation of the chemi- 
cal potential, 



N 



£ 1 



1 



2(nf-nD 



e£[A 2 + 4(e£ + e|)] (eg + e*)(fi^ + V 2 ) 



e£[A 2 + e|(4 + 4)] 



.(26) 



For fixed chemical potential we can find A(= A aa ) ver- 
sus l/fci?a s for different values of ^i, as shown in Fig. [5] 
For fixed number of particles equations [2U and must 
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be solved self-consistently. Fig. [7] shows the behavior 
of the superconductor order parameter and the chemical 
potential as function of the coupling parameter at T = 0. 
One can sc that in the BCS limit the chemical potential 
practically does not differ from the Fermi energy, and the 
superconducting gap is much smaller than Ep. Increas- 
ing the coupling, the pairs become more tightly bound 
and the chemical potential decreases and becomes nega- 
tive. 




-| . 1 . 1 . 1 1 1 . 1 r-^T " ^ 1 1 — "i 1 1 r 

-0.88 -0.86 -0.84 -0.82 -0.80 -0.78 -0.76 -0.74 -0.72 -0.70 -0.68 -0.66 
1/k a 



FIG. 6: (Color online) Zero temperature gap A as function 
of the interaction parameter l/kpa 3 , for the mass ratio a = 5 
and hybridization V=0.01 for different values of fixed chemical 
potential. 



FiglU The critical value of hybridization increases with 
the coupling parameter and tends to a saturation value. 
Notice that in this case there is a critical value of hy- 
bridization only for strong couplings as the chemical po- 
tential becomes negative. It can be understood analyzing 
Fig. O in which we find that the order parameter goes 
to zero for a definite value of l/fepOs just for negative 
values of the chemical potential (solid lines). 




0.6 n 



1/k a 



FIG. 8: (Color online) Critical values of the hybridization V c 
and chemical potential at T — as functions of the coupling 
strength l/kpa s , for the mass ratio a = 5. 




FIG. 7: (Color online) Zero temperature gap A and chemical 
potential fi as functions of the interaction parameter l/fcj?a s , 
for the mass ratio a — 5 and hybridization V=0.01. 



The behavior of the critical value of the hybridization 
as function of 1/kFds is shown in Fig. HI An increase 
in hybridization for a fixed value of the interaction sup- 
presses superconductivity continuously, as can be seen in 



UEMURA RELATION AND SCALING CLOSE 
TO THE SQCP 

Close to a superconducting quantum critical point, the 
vanishing of the superfluid density at T — is governed 
by the scaling relation [l7|, ps(T = 0) oc \g\ 2 ~ a ~ 2l/ where 
a and v are critical exponents associated with the SQCP 
and \g\ — \V — Vc\ measures the distance to the transi- 
tion. 

We are interested in dimensions d > 2, such that, for 
the superconducting quantum phase transition consid- 
ered here which is in the universality class of the density- 
driven Bose-Einstein condensation with dynamic expo- 
nent z = 2, the effective dimension d e ff = d + z >= 4. 
This coincides or is larger than the upper critical dimen- 
sion d c = 4 for this transition |I5]. In this case using 
the mean-field values, a = and v = 1/2, we find 
that ps scales linearly with the distance to the SQCP, 
p s (r = 0)cx \g\ as 3^0. 

On the other hand close to the SQCP, as \g\ — > 0, 
the superconducting critical temperature vanishes with 
the distance to the SQCP as, T c oc \g\v where ip is the 
shift exponent. For d + z > d c = 4, the interaction u 
between the bosons is dangerously irrelevant and the shift 
exponent is determined by the dynamic exponent and the 
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dimension of the system 15|, |l7|, |35j, i/; = z/(d + z — 2) 

a n 

Finally, we can derive a Uemura relation [36| relating 
the T — superfluid density and the critical temperature 
for close to a SQCP. Due to the proximity to the SQCP 
this is completely determined by scaling. For d + z > d c 
this is given by, T c oc pg with ip = z/(d + z — 2). In our 

2/3 

case with z = 2, we find, T c oc ps in <i = 2 and T c oc p s ' 
in d = 3. 



CONCLUSIONS 

We have studied the normal-superconductor quantum 
phase transition induced by hybridization in a two-band 
system using a new method recently introduced to deal 
with systems which are coupled to a space and time de- 
pendent perturbation. The results obtained are valid to 
first order in the perturbation and coincide with those of 
linear response theory. However, in our approach the ba- 
sic elements to be calculated are single particle Green's 
functions and not the usual two particle propagators of 
linear response theory. In the case of superconductivity, 
to calculate the superconducting response, we have to in- 
troduce a fictitious field which couples to the supercon- 
ductor order parameter. Starting from the normal phase, 
we obtained a generalized wave- vector and frequency de- 
pendent susceptibility. In the static and homogeneous 
limit, we argued that the divergence of this susceptibility 
signals the instability of the system to a superconduct- 
ing state. For the standard BCS problem with a single 
attractive band, this condition is similar to the Thou- 
less criterion for BCS superconductivity. For multi-band 
superconductors with intra-band interactions, we had ob- 
tained previously using a mean-field BCS approximation 
that superconductivity can be suppressed by increasing 
hybridization . Here we further clarify the mechanism 
of suppression and extend the mean-field approach to in- 
clude fluctuations close to the quantum superconductor- 
normal phase transition. We find that for finite attractive 
interactions, as hybridization increases, the BCS super- 
conductor has a smooth crossover to a superconducting 
breached pair-like state before becoming normal. The 
quantum superconductor-normal phase transition is in 
the universality class of the density-driven Bose-Einstein 
condensation. 

We have also investigated the case of fixed number of 
particles in the presence of hybridization. We have de- 
termined the T = phase diagram of the critical value 
of the hybridization versus the strength of the interac- 
tion given in terms of the scattering length. In this case 
superconductivity is also suppressed with the increase of 
hybridization, but only in the strong coupling limit for 
which the chemical potential becomes negative. 
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